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I  Introduction 


It  is  well  known  that  ray  optical  solutions  fail  at  and  near  caustics  of  ray  sys¬ 
tems  thus  requiring  one  to  modify  the  ray  solutions  in  such  regions  using  uniform 
asymptotic  techniques.  In  this  paper,  uniform  asymptotic  high  frequency  field  repre¬ 
sentations  are  presented  which  remain  valid  within  the  transition  regions  adjacent 
to  a  smooth  caustic  formed  by  the  envelope  of  the  geometrical  optics  (GO)  rays 
reflected  from  two  and  three-dimensional  smoothly  indented  perfectly-conducting 
boundaries  in  free  space.  These  field  representations  are  uniform  in  the  sense  that 
they  are  bounded  and  continuous  at  the  caustic;  furthermore,  outside  the  caustic 
transition  regions  they  recover  the  GO  reflected  ray  fields  which  are  valid  on  the  lit 
side  of  smooth  caustics,  and  they  properly  reduce  to  the  complex  ray  fields  given 
recently  by  Ikuno  and  Felsen  (1|  on  the  dark  side  of  such  caustics  wherein  real  rays 
do  not  penetrate. 

While  the  subject  of  evaluating  fields  in  caustic  regions  is  not  new,  the  uniform 
solutions  generally  presented  in  the  recent  past  mainly  emphasize  the  mathematical 
aspects  of  the  problem;  hence,  one  of  the  primary  goals  of  this  paper  is  to  try  to 
present  the  uniform  results  obtained  here  in  a  fashion  that  would  be  particularly 
useful  for  engineering  applications  requiring  a  ray  analysis  of  the  problems  of  elec¬ 
tromagnetic  (EM)  scattering  by  boundaries  that  can  generate  smooth  caustics  of 
reflected  rays.  Consequently,  an  important  feature  of  the  present  results  is  that  they 
provide  a  relatively  simple  prescription  not  only  for  evaluating  the  fields  on  the  lit 
side  of  smooth  caustics,  but  also  for  evaluating  the  fields  on  the  dark  side  of  such 
caustics.  It  is  noted  that  most  other  uniform  solutions  which  recover  the  GO  fields 
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on  the  deep  lit  side  of  the  caustic  appear  to  be  more  complicated  for  evaluating 
the  fields  within  and  outside  the  transition  region  on  the  dark  side  of  the  caustic; 
this  is  because  the  parameters  present  in  thr)se  solutions  typically  must  be  found 
from  a  knowledge  of  the  radius  of  curvature  of  the  caustic  and  the  perpendicular 
distance  from  the  caustic  to  the  observation  point.  While  it  is  known  that  it  also 
possible  to  more  directly  find  the  parameters  in  those  solutions  for  the  lit  side  from 
just  a  knowledge  of  the  GO  fields  which  need  to  be  corrected  within  the  caustic 
regions,  no  such  direct  simplification  appears  possible  in  these  solutions  to  find  the 
parameters  for  evaluating  the  fields  on  for  the  dark  side  of  the  caustic.  On  the  other 
hand,  the  parameters  occurring  in  the  present  solutions  can  be  found  relatively  eas¬ 
ily  from  the  stationary  phase  points  in  the  physical  optics  (PO)  approximation  for 
the  radiation  integral  over  the  currents  induced  on  the  reflecting  boundary  by  the 
incident  field.  It  is  noted  that  the  PO  approximation  employs  the  usual  geometrical 
optics  (GO)  field  to  find  these  induced  currents.  This  PO  approximation  is  valid 
whenever  the  smooth  reflecting  boundary  is  electrically  large  and  well  illuminated; 
these  conditions  are  assumed  to  hold  in  the  present  analysis  and  are  the  same  as  the 
ones  required  in  the  GO  reflection  calculations.  There  are  two  real  stationary  phase 
points  in  the  PO  radiation  integral  over  the  reflecting  boundary  for  an  observation 
point  on  the  lit  side  of  a  smooth  caustic  of  the  reflected  rays;  these  stationary  points 
coalesce  to  form  a  second  order  stationary  point  for  observation  points  on  the  caus¬ 
tic  itself,  and  then  they  become  two  complex  stationary  points  (which  are  complex 
conjugates  of  each  other)  for  observation  points  on  the  dark  side  of  the  caustic  even 
though  only  one  of  these  eventually  contributes  to  the  fields  on  the  dark  side  outside 
the  caustic  transition  region.  As  usual,  the  relevant  stationary  phase  points  can  be 
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found  simply  by  requiring  the  spatial  derivatives  of  the  phase  function  in  the  PO 
radiation  integral  over  the  reflecting  surface  to  vanish  as  described  later  in  Section 
II. 

It  is  useful  at  this  juncture  to  briefly  review  some  of  the  previous  related  work. 
Rahnavard  and  Rusch  [2]  treated  the  problem  of  electromagnetic  (EM)  scattering 
by  a  shaped  subreflector  surface  with  a  point  of  inflection.  Their  PO  based  asymp¬ 
totic  solution  is  valid  in  the  transition  region  adjacent  to  the  locally  smooth  portion 
of  the  reflected  ray  caustic  which  results  from  the  effect  of  an  inflection  point  on 
the  subreflector  geometry.  However,  their  solution  is  only  partially  uniform  as  it 
does  not  recover  the  GO  reflected  field  on  the  lit  side  of  the  caustic;  also,  they 
introduced  an  additional  approximation  in  the  asymptotic  treatment  of  the  PO  ra¬ 
diation  integral  for  evaluating  the  field  on  the  dark  side  of  the  caustic  instead  of 
making  use  of  the  complex  stationary  points  in  the  integrand.  Nevertheless,  their 
solution  provided  some  useful  numerical  results  for  shaped  subreflector  design  when 
it  was  employed  within  its  range  of  validity.  Albertsen,  et  al.  [3]  obtained  a  solution 
valid  near  smooth  caustics  of  edge  diffracted  rays;  they  performed  an  asymptotic 
evaluation  of  the  radiation  integrals  over  equivalent  edge  currents  to  arrive  at  their 
solution.  Next,  they  introduced  a  heuristic  factor  in  their  solution  so  that  the  result¬ 
ing  modified  solution  would  reduce  uniformly  to  the  ray  solution  (based  on  Keller’s 
edge  diffracted  ray  fields)  on  the  lit  side  of  the  caustic;  they  also  introduced  the 
same  approximation  as  in  [2]  for  the  dark  side  of  the  caustic.  Kravatsov  [4],  and 
Ludwig  [.5|,  independently  showed  that  the  leading  terms  of  an  asymptotic  solution 
which  remained  uniformly  valid  across  a  smooth  caustic  contained  not  just  an  Airy 
function  as  in  [2,3],  but  also  its  derivative  which  was  missing  in  [2,3].  Ludwig’s 
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ansatz  [5]  is  actually  based  on  a  uniform  asymptotic  evaluation  of  an  integral  with 
two  nearby  saddle  points  as  given  by  Cdiester,  Friedman  and  Ursell  [6].  The  lat¬ 
ter  asymptotic  evaluation  is  summarized  lucidly  in  texts  by  Felsen  and  Marcuvitz 
[7],  and  Stamnes  [8].  While  the  term  involving  the  derivative  of  the  Airy  function 
vanishes  at  the  caustic,  it  is  of  the  same  order  as  the  Airy  function  term  exterior 
to  the  caustic  transition  region  and  its  inclusion  is  therefore  necessary  to  uniformly 
recover  the  corresponding  ray  solution  deep  in  the  lit  region.  It  is  noted  that  the 
solutions  developed  in  this  paper  also  contain  the  Airy  function  and  its  derivative 
since  the  integrals  leading  to  these  solutions  are  evaluated  asymptotically  using  the 
techniques  developed  by  Chester,  et  al.  [6,7].  Ludwig’s  analysis  [5]  provides  im¬ 
portant  extensions  to  the  earlier  work  of  Kay  and  Keller  [9],  and  in  particular  of 
Buchal  and  Keller  [lOj.  Ludwig’s  approach  [5]  requires  one  to  initially  assume  a 
general  form  of  the  uniform  solution  given  by  Chester,  et  al.  [6j  and  the  amplitudes 
and  arguments  of  the  uniform  functions  (involving  the  Airy  function  and  its  deriva¬ 
tives)  are  then  obtained  by  requiring  the  solution  to  satisfy  the  wave  equation  in  the 
high  frequency  limit.  However,  this  systematic  and  elegant  approach  for  the  two- 
dimensional  case  by  Ludwig  [.5]  appears  to  be  cumbersome  especially  for  treating 
the  corresponding  three-dimensional  ca.se.  More  recently,  phase-space  techniques 
have  been  employed  by  Maslov  [11,12]  for  evaluating  the  fields  in  the  neighborhood 
of  ray  caustics;  in  that  method,  the  field  is  eventually  expressed  as  a  spectral  inte¬ 
gral  which  can  be  transformed  locally  into  a  canonical  form,  as,  for  example,  the 
one  treated  by  Chester,  et  al.  [6]  for  the  ca.se  of  a  smooth  caustic.  The  general  re¬ 
lationships  between  caustics,  their  as.sociated  field  singularities  and  ray  geometries 
are  also  expressed  in  a  systematic  fashion  in  catastrophe  theory  [13],  which  defines 
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caustics  as  singularities  of  ray  families  and  classifies  caustics  according  to  their  sin¬ 
gularities.  In  particular,  the  smooth  caustic  is  classified  as  a  “fold”  catastrophe. 
Catastrophe  theory  provides  a  mapping  for  the  phase  function  in  an  integral  repre¬ 
sentation  for  the  field  near  a  caustic  so  that  it  can  like  wise  be  locally  transformed 
into  a  canonical  integral.  Very  recently,  Arnold  [14,15]  has  essentially  extended 
Maslov’s  method  to  include  discontinuities  in  a  scattering  boundary;  there,  he  em¬ 
ploys  a  spectral  synthesis  of  uniform  wavefunctions  for  analyzing  fields  near  smooth 
caustics  based  on  ideas  common  to  Maslov’s  method  and  catastrophe  theory.  As 
in  the  present  approach,  Arnold’s  method  is  readily  applicable  to  treat  the  case  of 
smooth  caustics  in  three  dimensions.  However,  even  though  these  previous  elegant 
treatments  in  [4]-[6],[10]-[15]  provide  the  general  concepts  involved  for  analyzing  the 
fields  in  caustic  regions,  they  are  presented  at  a  relatively  sophisticated  mathemat¬ 
ical  level  and  the  final  solutions  which  generally  avoid  any  explicit  treatments  of 
the  fields  on  the  dark  side  of  the  caustic  do  not  appear  to  be  readily  ammenable  for 
use  in  engineering  applications.  These  treatments  in  [4]-[6],[10]-[15]  provide  uniform 
solutions  (usually  explicitly  for  the  lit  side  of  the  caustic)  which  contain  parameters 
that  are  essentially  expressed  in  terms  of  the  caustic  (or  corresponding  wavefront) 
geometry;  hence,  it  is  necessary  first  to  know  the  details  of  the  caustic  geometry 
especially  for  evaluating  the  fields  on  the  dark  side  of  the  caustic.  The  latter  com¬ 
plication  is  avoided  in  the  present  approach.  It  is  noted  that  the  present  approach 
can  be  easily  generalized  to  deal  with  the  reflection  from  an  impedance  boundary 
or  from  a  dielectric  interface. 
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The  problems  of  two  and  three-dimensional  smooth  caustics  are  analyzed  in 
essentially  the  same  fashion  in  the  present  approach.  If  there  are  two  caustics,  the 
present  analysis  is  valid  across  each  caustic  provided  that  the  obervation  point  is 
far  from  the  region  where  the  two  caustics  merge  together  to  from  a  cusp.  Near 
the  cusp,  a  different  approximation  is  necessary.  There  are  generally  two  smooth 
caustics  present  in  the  (3-D)  case;  the  present  analysis  is  valid  across  each  smooth 
caustic  for  the  latter  case  provided  that  these  two  caustic  surfaces  do  not  intersect. 
In  practical  applications,  the  surfaces  which  scatter  are  of  finite  extent;  hence, 
there  may  be  situations  where  the  caustic  can  terminate  in  which  case  the  present 
analysis  is  valid  across  a  locally  smooth  portion  of  the  caustic  which  is  sufficiently 
far  from  the  caustic  termination,  and  ‘he  g"*nerally  smaller  effect  of  that  termination 
must  be  found  separately.  The  uniform  results  are  developed  and  presented  here 
in  Sections  II  and  III  for  the  2-D  and  3-D  vector  EM  cases,  respectively,  with 
all  the  parameters  being  defined  explicitly  for  both  the  lit  and  the  dark  side  of 
the  caustics  in  a  form  suitable  for  engineering  applications.  These  uniform  results 
are  employed  in  this  paper  to  analyze  the  scattering  by  boundaries  which  contain 
inflection  points.  In  particular,  numerical  results  have  been  obtained,  as  shown  in 
Section  IV,  for  the  far  zone  EM  plane  wave  backscattering  by  a  concave-convex 
shaped  boundary  with  an  edge,  and  also  by  a  smoothly  indented  cavity,  both  of 
which  contain  points  of  inflection  on  their  boundaries.  The  numerical  results  for  the 
concave-convex  geometry  with  an  edge  are  shown  to  compare  well  with  independent 
moment-method  calculations.  Numerical  results  for  the  smoothly  indented  cavity 
are  compared  with  those  of  a  rectangular  cavity  of  the  same  length  and  depth.  It  is 
seen  that  the  scattering  from  the  smoothly  indented  cavity  can  be  of  the  same  order 
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within  the  caustic  regions  as  that  from  a  non-smooth  rectangular  shaped  cavity. 
An  time  dependence  is  assumed  and  suppressed  in  the  following  analysis. 


II  An  Asymptotic  Solution  for  the  Field  Associ¬ 
ated  with  a  Two  Dimensional  Smooth  Caustic 


An  asymptotic  high  frequency  analysis  is  presented  in  this  section  for  predicting 
the  EM  fields  within  the  caustic  regions  of  rays  reflected  from  smoothly  indented 
perfectly-conducting  two  dimensional  boundaries  in  free  space.  The  excitation  in 
these  problems  is  assumed  to  be  produced  by  a  localized  source  which  is  sufficiently 
far  (in  terms  of  the  electrical  wavelength)  from  the  boundaries  so  that  the  field  inci¬ 
dent  from  the  source  can  be  represented  ray-optically.  If  the  excitation  is  due  to  an 
electrically  large  or  extended  source  distribution,  then  one  can  quantize  that  source 
distribution  so  that  each  quantized  source  element  satisfies  the  above  criterion;  the 
complete  reflected  field  in  this  case  can  be  found  by  a  superposition  of  the  fields 
reflected  from  the  boundaries  when  they  are  illuminated  by  each  of  the  quantized 
source  elements. 

The  scattering  by  geometries  depicted  in  Figures  1(a)  and  (b)  represent  typical 
situations  of  interest  here  in  which  smooth  caustics  of  reflected  rays  are  formed. 
Both  the  non-uniform  as  well  as  the  uniform  asymptotic  high  frequency  approxi¬ 
mations  of  the  physical  optics  integral  for  the  scattered  field  in  two-dimensions  are 
presented  separately  in  Parts  A  and  B  of  this  section,  respectively.  The  non-uniform 
asymptotics  furnish  the  usual  pair  of  GO  reflected  ray  fields  on  the  lit  side  of  the 
smooth  caustic  of  reflected  rays;  these  two  real  reflected  rays  are  associated  with  the 
two  real  stationary  phase  points  in  the  PO  integrand,  and  they  become  degenerate 
as  the  stationary  points  coalesce  for  observation  points  on  the  caustic  where  their 
associated  field  exhibits  a  singularity.  The  non-uniform  asymptotics  also  provides 
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a  result  for  the  field  on  the  dark  side  of  the  caustic  which  is  associated  with  only 
one  of  the  two  complex  stationary  points;  as  indicated  earlier  by  Ikuno  and  Fel- 
son  [1],  this  evanescent  (exponentially  decaying)  field  can  be  described  in  terms  of 
a  complex  ray  specularly  reflected  from  a  complex  extension  of  the  surface.  The 
other  complex  (conjugate)  stationary  point  gives  an  unbounded  result  and  it  must 
be  discarded  in  the  non-uniform  analysis  on  physical  grounds  as  its  contribution 
violates  the  radiation  condition.  The  non-uniform  result  for  the  dark  side  also  ex¬ 
hibits  a  field  singularity  at  the  caustic  as  does  the  GO  result  for  the  lit  side.  It  is  for 
the  sake  of  completeness  that  the  non-uniform  asymptotic  treatment  is  included; 
moreover,  it  serves  to  introduce  the  notation  which  is  also  used  later  in  the  uniform 
asymptotic  analysis.  Of  course,  the  non-uniform  result  is  not  valid  at  and  near  the 
caustic;  instead,  one  must  employ  a  uniform  solution  there  as  discussed  in  Part  B. 

II-A  Non-Uniform  2-D  Analysis: 

The  electric  field  E’{P)  at  any  point  P  exterior  to  the  boundary  C  of  Figures  1(a) 
and  1(b)  which  here  constitutes  the  field  scattered  by  the  boundary  can  be  expressed 
in  terms  of  the  usual  radiation  integral  over  the  electric  current  J  induced  on  C  by 
the  source  at  P'  as 

IcZ  r 

E‘{P)  -  /  [ff  '<  ff  X  JiQ')]  Hi^\kR)de[  (1) 

where  J(Q')  is  the  value  of  J  at  any  point  Q'  on  the  boundary  C  and  R  =  Q'P. 
Also,  the  point  P  cannot  be  too  clo.se  to  the  reflecting  boundary  for  (1)  to  be  valid. 
Also  k  and  Z„  are  the  wavenumber  and  impedance  of  free  space  respectively  and 
H]^\kR)  is  the  usual  cylindrical  Hankel  function  of  the  second  kind  of  order  zero. 
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To  be  specific,  P  =  Pi  for  an  observation  point  on  the  lit  side  of  the  caustic  and 
P  —  Pg  for  an  observation  point  on  the  dark  side  of  the  caustic,  respectively.  In 
the  PO  approximation  to  (1),  the  J{Q')  is  assumed  to  be  given  by  GO  as  follows; 


J{Q')  « 


2h  X  H'{Q'),  on  the  lit  portion  of  the  boundaries 
0,  on  the  shadowed  portion  of  the  boundaries 

(where  the  direct  illumination  from  the  source  is  absent) 


(2) 


in  which  H'{Q')  is  the  magnetic  field  at  Q'  which  is  directly  incident  from  the 
source  whose  phase  center  is  at  P',  and  n  is  the  unit  normal  vector  to  the  boundary 
C  at  {Q')  as  shown  in  Figures  1(a)  and  1(b).  Under  the  previously  mentioned 
assumptions,  the  H'{Q')  may  be  expressed  ray  optically  as 


ff'(Q') 


z:^ X 


;  with  i'o  -.4  =  0 


(3) 


It  is  also  understood  that  H*{Q')  ~  Z~^s*  x  E'{Q'),  in  which  s'  =  P'Q',  s'  =  ||fj, 
and  E'{Q')  is  the  electric  field  associated  with  H'(Q')\  it  is  clear  from  (3)  that 
E'{Q')  =  A^ .  Using  the  large  argument  approximation  of  H^\kR)  for  kR  >> 
1,  and  the  assumed  current  in  (2)  within  (1)  yields  the  PO  integral  for  E‘{P)  at  an 
observation  point  P  as 


E-(P) 


Ijk 

2ir 


[R  X  R  X  {n  X  s'  X  ^)] 


-jk(s'+R) 


vTr 


(4) 


The  integral  in  (4)  can  be  evaluated  for  large  k  via  the  method  of  stationary  phase 
[7,8].  The  stationary  phase  points  are  determined  by  the  condition 


i[  ■  V[s'  +  R)^0 


(5) 
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ill  which  along  the  two-dimensional  reflecting  boundary  C.  Both  s'  and 

R  are  functions  of  t'  (arc  length  to  Q'  along  C  from  some  reference  point  O  on  C). 
The  solution  to  (5)  yields  the  stationary  phase  points  at: 


and  ,  at  Qa  and  Qt,  respectively  on  C, 
when  P  =  Pi  for  an  observation  point  on  the  lit  side  of  the  caustic 


(6) 


where  for  the  sake  of  being  specific  one  requires  the  following,  in  addition  to  (5): 


d^(s'  +  R) 


di'^ 


>0; 


+  R) 


-^la 


<0 , 


d-d 


and 


and  at  Ql  and  respectively, 
when  P  =  Pi  for  an  observation  point  on  the  dark  side  of  the  caustic. 


The  points  Q'  and  Ql  are  associated  with  the  complex  stationary  phase  points 
and  (-1  which  are  conjugates  of  each  other;  thus,  and  Ql  must  lie  on  the  complex 
extension  of  the  boundary  C  as  indicated  in  [1].  At  the  stationary  points  which 
correspond  to  points  of  specular  reflection  at  Qa  and  Qb  on  C,  for  P  =  Pi,  the 
following  conditions  hold;  namely: 


R  X  n  =  s' 


X 


n 


(8) 


R  ■  fi  =  —  s'  •  n  =  cos  9' 


(9) 


JR  —  s'  =  2(n  •  R)h  —  2  cos  9'h 


(10) 
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for  Q'  =  b 


(11) 


for  Q'  =  Q^  b 


h  =  fia  i,  for  Q'  =  Q^  b 


O'  = 


for  Q'  =  Qa,b 


In  the  neighborhood  of  the  stationary  points  and  Qb  when  P  =  Pi,  one  arrives 
at  the  following  local  approximation  for  the  phase  term  k{s'  +  R)  as: 


s' +  R^  sib  +  ^,6  +  TyCMQ..^) cos^olbit,  - e[„,,bf 


for  Q'  near  Qa,b  • 


(12) 


At  high  frequencies;  i.e.,  for  large  k,  the  dominant  contribution  to  (4)  comes  from 
the  stationary  phase  points  Qa  and  Qb  on  C  when  P  =  Pi;  thus,  (4)  may  be 
approximated  via  (12)  as  usual  to  obtain: 


E’{Pi) 


X  X  (n„  X  i;  X  A{Q„))]  ^ 


Jca 

+  1^6  X  X  (nf,  X  i;  X  i4(Qi,)))' 

.jT  df;]| 


\/^'b^b 


(13) 
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in  which  A{Qa,b)  denotes  the  value  of  A  at  Qa,b‘,  also  Gi{Qa,b)  is  defined  as 

Gi{Qa,b)  =  —  i — ;  ,  at  Qa,b  • 

^a,b  Pla.lb 


(14) 


It  can  be  shown  that 


1 


+ 


(15) 


P'la.lb  <,6  PgiQafi)  COS  dlj, 

The  domains  of  integration  and  Cb  include  the  isolated  or  disjoint  neighborhoods 
of  the  points  Qa  and  Qb  on  C.  The  Pg(Qa,b)  denotes  the  radius  of  curvature  of 
the  boundary  C  along  at  Qa,bi  furthermore,  corresponds  to  the  radius  of 

curvature  of  the  reflected  wavefront  associated  with  the  GO  rays  reflected  from  the 
points  Qa,b  on  C,  and  j,  is  the  usual  angle  of  incidence  (made  by  the  incident  ray 
and  the  normal  to  C  at  Qa_b)-  An  evaluation  of  the  integrals  over  Co  and  Cb  in  (13) 
after  extending  the  limits  of  integration  to  infinity  via  the  stationary  phase  approach 
(which  yields  ^  wherein  5(gin(Gi(Qo))  =  +1  and 

Sgn(Gi(Qb))  =  —1  for  Gi(Qa)  >  0  and  Ci(Qfr)  <  0,  respectively.  )  leads  to 


E^(Pl) 


Ijk 

2ir 


[R^  X  R„x  (h„  X  j;  X  A(Qo))] 


27r 


kGiiQ a)  COS^ 

+  [^i,  X  Rb  X  {hb  X  Sfc  X  i4(Qe.))] 


^-j’,Sgn{G(Qa)) 


I  2it 

kGi(Qb)cos'‘  $1 


\f^b 

^-jfS»n(G(Qk))| 


(16) 


It  is  noted  that 


Ra^b  ^  Ra,b  ^  i.'^a,b  ^  ^a,6  ^  A(^o,6))  —  ( ^"^a.b^a^b  '  A(^o,6)  ^OS  ^o,6  (10 
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in  which  I  is  an  identity  dyad  (so  that  I  ■  A  —  A).  In  two-dimensions,  where  one  is 
restricted  to  the  z  =  0  (say)  plane;  i.e.,  to  the  (x,y)  plane,  then  /  =  xi  -f-  yy.  For 
example,  /  =  xx  -f-  yy  -|-  zz  in  a  rectangular  coordinate  system.  Thus,  incorporating 
(14)  and  (17)  into  (16)  yields: 


(18) 


where 


El.biPl)  =  (2ha.6na,(,  -  i)  •  A{Qa,b) 


Pu.lb 


\  Pla.Xb  +  '*a,fc 


(19) 


Since,  A(Qa.b)^  7T=  identified  as  E’{Qa_i,}  (see  comments  below  (3)),  the 


preceeding  result  becomes: 

EI,{Pl)  =  E\Q.,b)-R{Qa.b). 

where 


Pla.lb 


\  Pu.Xb  +  «a,(, 


(20) 


R(Qa.b)  ^‘^na.bfla.b  -  I  (21) 

can  be  recognized  as  the  dyadic  reflection  coefficient  for  a  perfectly  conducting 
boundary.  Consequently,  the  non-uniform  asymptotic  result  for  E‘{Pl)  in  (18) 
together  with  (20)  can  be  identified  as  the  sum  of  the  fields  associated  with  the 
two  GO  rays  reflected  from  the  points  Qaj,  as  shown  in  Figures  1(a)  and  1(b).  It 
is  noted  that  one  must  require  Gi{Qa)  >  0  and  Gi{Qi,)  <  0  to  be  consistent  with 
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Figures  1(a)  and  1(b);  thus, 


and 


since  G'i(Qa)  >  0  , 


(22) 


e-'*  ,  since  Gi(Q6)  <  0  ,  (23) 

must  be  true  in  (19)  and  (20)  as  required  by  (16).  Clearly,  (22)  and  (23)  indicate 
that  the  principal  (positive)  branch  of  the  square  roots  must  be  taken.  As  usual, 
the  additional  phase  e-'  j  in  (23)  indicates  that  the  GO  ray  reflected  from  the  point 
Qk  touches  the  caustic  before  reaching  Pi  as  shown  in  Figures  1(a)  and  1(b).  It  is 
convenient  to  express  the  reflection  coefficient  R  in  terms  of  the  unit  vectors  fixed 
in  the  incident  and  reflected  rays  [16]  as  shown  in  Figure  1(c)  ;  thus. 


R  —  RtCiii  +  /Zfceljejj  ; 

iZj  =  ^1  for  a  conducting  boundary.  (24) 

h 

The  non-uniform  asymptotic  high  frequency  or  the  GO  approximation  for  the 
refelcted  field  E’(Pl)  at  Pi  which  is  valid  on  the  lit  side  of  the  caustic  but  external 
to  the  caustic  transition  region  is  given  explicitly  via  (16)  through  (24)  as: 

E‘(Pi)^  E\Q,)  -Mq^) 

^E'(Q,)-R{Qi,)-j\J-7^  •  (25) 

IV  + 

All  the  quantities  in  (25)  have  been  defined  previously.  It  is  noted  that  R{Qa,b) 
in  (25)  is  as  given  in  (24)  except  that  the  unit  vectors  cx>  c||  evaluated 


15 


separately  for  the  first  and  second  terms  on  the  right  hand  side  of  (25)  corresponding 
to  each  of  the  rays  reflected  from  Qa  and  Q(,,  respectively. 

At  the  caustic  aj  =  pj,  and  =  —  pL  !  hence  the  GO  reflected  field  E*{Pl) 

b  lb  lb  lb 

in  (25)  becomes  singular  there.  The  representation  in  (25)  is  not  valid  at  and  near 
such  a  caustic  of  reflected  rays. 

One  can  likewise  extend  the  stationary  phase  analysis  into  the  shadow  region; 
in  the  latter  case,  one  can  show  that  the  scattered  field  E^{Ps)  at  a  point  Ps  in  the 
shadow  region  is  given  by: 

E‘{P.)  ~  E'(Ql)  ■  l{Qt)J  /f  ,,  (26a) 

y  Pb  +  ■*b 

with 

E'iQl)  =  (26b) 

The  contribution  in  (26)  is  associated  with  the  complex  stationary  phase  point 
QJ  corresponding  to  i'  =  t'f,  (see  (7)),  which  must  lie  on  the  complex  extension  of 
the  reflecting  boundary  as  discussed  by  Ikuno  and  Felsen  [2];  the  contribution  from 
the  other  (conjugate)  stationary  phase  point  at  Q'  corresponding  to  C  =  must 
be  discarded  on  physical  grounds  since  it  gives  rise  to  an  exponentially  increasing 
field  amplitude  as  mentioned  before.  The  contribution  in  (25)  can  be  shown  to  be 
exponentially  decreasing  corresponding  to  an  evanescent  field  at  Ps  on  the  dark  side 
of  the  caustic.  The  superscript  c  in  the  quantities  appearing  in  (26)  denotes  that 
those  quantities  are  associated  with  the  “complex”  stationary  point  Ql- 

It  is  noted  that  the  quantity  pj'  on  the  right  hand  side  of  (26)  is  calculated  via 
an  analytic  continuation  of  the  corresponding  expression  for  pj  in  (15).  This  follows 
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because  the  point  of  specular  reflection  on  the  actual  boundary  in  real  space  for 
an  observation  point  P  —  Pl  on  the  lit  side  of  the  caustic  is  analytically  continued 
to  a  complex  location  which  obviously  must  lie  on  the  complex  extension  of 
the  original  boundary  when  P  =  Ps  for  an  observation  point  on  the  dark  side  of 
the  caustic  [1].  Thus,  when  the  right  hand  side  (RHS)  of  the  expression  in  (15)  is 
evaluated  at  the  real  stationary  point  Q(„  then  one  obtains  the  value  denoted  by 
pI  which  appears  on  the  left  hand  side  (LHS)  of  (15);  likewise,  when  the  (RHS)  of 
(15)  is  evaluated  at  a  complex  stationary  point  then  one  obtains  the  analytic 
continuation  of  pj  from  its  value  at  the  real  coordinates  of  to  its  value  at  the 
complex  coordinates  of  Ql-  The  manner  in  which  the  quantities  s",  Pg{Ql), 
cosd'i^  (=  value  of  cos^*  in  (10)  evaluated  at  Q^),  etc.,  are  evaluated  for  the 
complex  stationary  point  by  a  simple  analytic  continuation  of  the  corresponding 
expressions  given  in  (8)-(ll)  and  (15)  for  real  stationary  points,  is  described  in  the 
Appendix  for  the  sake  of  completeness  and  convenience.  Also,  the  H(Q^)  may  be 
found  by  analytic  continuation  of  R(Qb)  in  (21)  to  the  complex  location  Ql  in  place 
of  Qb,  or  by  evaluating  R  in  (24)  at  Ql;  either  of  these  procedures  should  yield  the 
same  result  for  i2(Qj).  If  (21)  is  used,  then  the  nj,  therein  which  is  evaluated  at  Qb 
must  be  replaced  by  hi  when  it  is  evaluated  at  similarly,  if  one  employs  (24),  then 
the  e\,  e||,  e|j  must  be  replaced  by  their  analytically  continued  values  e\,  eJI"  and  elj' 
for  the  complex  stationary  point  as  shown  in  the  Appendix.  The  expression  in 
(26)  also  becomes  singular  at  the  caustic  where  Ql,  Pi^, s",  etc.,  become  real  and 
p^^  =  —  I  pj^  I  as  well  as  =|  p^^  |;  hence,  (26)  is  also  not  valid  at  and  near  the 
caustic.  A  uniform  representation  valid  at  and  near  the  reflected  ray  caustic  which 
recovers  the  results  in  (25)  and  (26)  for  the  lit  and  dark  sides  of  the  caustics  that 
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lie  outside  the  caustic  transition  region  is  described  next. 

II- B  Uniform  2-D  Analysis 

The  starting  point  for  the  analysis  is  again  the  integral  representation  for  the  scat¬ 
tered  field  E‘{P)  which  is  given  in  (1).  This  integral  will  be  evaluated  via  a  uniform 
asymptotic  procedure  [6]  to  overcome  the  limitations  (near  the  caustic)  of  the  non- 
uniform  solution  which  was  developed  in  Part  A  of  this  section.  Again,  one  may 
make  the  usual  approximations  which  lead  to  (4)  as  before.  The  results  in  (5),  (6) 
and  (7)  are  still  valid.  However,  instead  of  incorporating  the  local  phase  approxima¬ 
tion  of  (12)  into  (4),  which  is  valid  if  the  stationary  points  are  isolated  (or  not  close 
together)  as  it  happens  when  one  is  outside  the  caustic  transition  region,  one  mast 
now  employ  a  different  procedure  which  remains  valid  even  where  the  stationary 
points  come  close  together  as  it  happens  for  observation  points  near  the  caustic. 
Thus,  instead  of  (12),  one  employs  the  following  transformation  when  P  =  Pi  on 
the  lit  side  of  the  caustic: 

3'  +  R=^-e-Cii  ;  O>0  (27) 

where 

^  iw +  »:)  +  (»;  m:)i  (28) 

and 

</  =  4  ((>:  +  »;)-(■< +  >:)i  •  (29) 
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all  the  quantities  on  the  RHS  of  (28)  and  (29)  have  been  defined  previously  in  Part 
A  (also  see  Figures  1(a)  and  1(b)).  Substituting  (27)  into  (4)  yields 

E‘(Pl)  «  r  dt  F(t)  (30) 

J -oo 

with 

[Tie  .  .  ^  ^  df 

(31) 

The  limits  of  integration  in  (30)  have  been  extended  to  infinity  even  if  C  is  a  finite 
closed  contour  because  in  the  present  development  one  is  interested  only  in  the  con¬ 
tribution  from  the  saddle  points  which  dominates  near  the  caustic;  the  contribution 
from  the  end  points  at  ±oo  is  vanishingly  small.  Any  other  contributions  to  the 
radiation  integral  in  (4)  which  may  arise  due  to  a  finite  closed  contour  G  must  be 
added  separately  and  are  not  dealt  with  here. 

Following  Chester,  et  al.  [6|,  one  may  now  expand  F{t)  in  (31)  as 

m  =  f:  [d^{p  -  Cer  +  -  or]  m 

m=0 

Retaining  only  the  leading  term  (m  =  0  term)  in  the  above  expansion  of  (32)  yields 
the  following  result  for  E’{Pi)  in  (30): 

^‘(Pl)  ~27rj  £-■'''<*  \Ptk  Ui{-kh()  i  j0'^Q(k-Ui'{-kHi)] 

+0{k-^)  (33) 

in  which  Ai{(r)  =  e^-’*  df  ;  and  .4?'(<t)  =  At(<r).  It  can  be  shown  that 

the  Pc  and  Qc  in  (33)  are  given  by: 

BaCos^O'^  * 

.  \/<^l(Qa)cOs2^j, 
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(34) 


BhCosOl 

\J4  y/C^iiQk)  cos  0i 

with 

Bo,(,COS^^  j,  =  Ra,b  X  fia.b  X  [ha,6  X  77^  5  X  ^(Qa.b)]  .  (35) 

It  follows  from  (17),  (19)  and  (20)  that  the  preceeding  result  in  (35)  can  be 
written  as 


Ba.b  cosd'^^ 


A(Q„,b)-R(Qa.b)  cos  91, 


(36) 


in  which  A{Qa,b)  's  given  by  E'{Qa,b)  =  AiQa.b)- 


-yfcj*  . 

■'  a,b 


as  indicated  previously  below 


(26).  Incorporating  (37)  together  with  (14)  into  (34)  and  (35)  gives 


p< 

U<J 

-  - 

2  V?r 

A(QJ-R(QJ  /  p\, 
S'  V  P\a  + 


e 


A{Qb)-RKQb)  rp^ 


\/»b 


P\b  +  K 


0+31 


(37) 


with  the  understanding  that  (22)  and  (23)  are  true  in  (37).  It  is  noted  that  v4(Qa_i,) 
can  depend  on  not  only  via  e||  and  ex  in  terms  of  which  the  polarization  of 
can  be  expressed,  but  there  could  be  additional  pattern  information  of  the  source 
contained  in  i4  as  a  result  of  which  A  can  also  depend  on  ,  (that  along  with  e||  and 
ex  in  turn  depends  on  Qa,b)-  On  the  lit  side  exterior  to  the  caustic  transition  region, 
the  result  for  E’(Pi)  in  (33)  toegther  with  (28),  (29)  and  (37),  uniformly  reduces  to 
the  GO  result  in  (25).  One  can  verify  this  GO  limit  by  employing  the  asymptotic 
approximations  for  Ai{~k3  (^()  and  Ai'{  in  (33)  for  G  >>  1  corresponding  to 
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the  lit  side  of  the  caustic  which  lies  exterior  to  the  caustic  transition  region;  these 
asymptotic  approximations  are; 

^  ^  smilkc}  +  ^)  ;  0  »  1  ,  (38) 

\/7r  U  4 

Ai'(-klCf)  - . 4"  kk;  cos(?k|  +  j)  ;  0  >>  1  .  (39) 

A  uniform  result  similar  to  that  in  (33)  can  also  be  obtained  for  the  shadow  side 
of  the  caustic;  this  uniform  analysis  for  the  shadow  side  is  presented  next.  Once 
again,  the  starting  point  is  the  integral  representation  in  (4).  A  transformation 
similar  to  that  in  (27)  is  introduced  into  the  phase  term  in  the  integrand  of  (4) 
when  P  =  Pi  on  the  shadow  side  of  the  caustic  as 

s*  =  J  ;  <,  >0,  (40) 

to  account  for  the  proximity  of  the  stationary  phase  points  and  Ql  in  the  phase 
function  (s'  4  /?),  which  occur  on  the  complex  extension  of  the  original  boundary 
(see  (7)).  Also, 

=  l  !(.< +  »;')  + (4 +  .C)I  HD 

and 

cl  =jI  [«  +  ^r)-(c  +  c)]  •  (^2) 

The  and  .sj''  in  (41)  and  (42)  are  the  same  as  those  seen  previously  in  (26a;  26b); 
these  quantities  are  the  analytically  continued  values  of  s),  and  from  those  asso¬ 
ciated  originally  with  the  real  stationary  point  Qi,  to  the  complex  stationary  point 
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at  Ql  as  shown  in  the  Appendix.  Likewise,  and  in  (41)  and  (42)  represent 
the  analytically  continued  values  of  and  from  those  associated  originally  with 
the  real  stationary  point  to  the  complex  stationary  point  at  Q'.  Substituting 
(40)  into  (4)  yields 

E’(Ps)  ^  r  dtF(t )€-■’''< (43) 
J  —  OO 

where  F{i)  was  defined  earlier  in  (31).  Again,  one  begins  by  expanding  F{t)  as  in 
(32);  next,  retaining  only  the  m=0  term  in  that  expansion  leads  to  the  following 
uniform  result  for  the  scattered  field  E‘‘{Ps)  valid  on  the  dark  side  of  the  caustic: 

E‘(Ps)--2irj  e-^'‘<‘[P,k-lAtiklC)  +  0"Q.k-lAi'(klC)] 

+Oik-^)  (44) 


in  which 


p. 

d  [k  \ 

I 

to 

a(q:)-h(q:) 


r  la 


MQI)-R{QI) 

\/4' 


p\l  + 


p\l 


p[b  + 


(45) 


It  is  noted  that  since  C  in  (42)  is  positive  (see  (40)),  the  quantity  within  the 
square  brackets  in  (42)  satisfies  the  conditions 


(46) 


/ml«  +  .,r  -  +  <')!  >  0 
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On  the  shadow  side  exterior  to  the  caustic  transition  region,  the  result  for  E‘{Ps) 
in  (44)  together  with  (41),  (42),  and  (45)  reduces  uniformly  to  the  analytically  con¬ 
tinued  GO  result  in  (26)  for  the  field  reflected  from  Ql-  One  can  verify  this  limiting 
value  in  (26)  for  the  expression  in  (44),  by  employing  the  asymptotic  approxima- 

3  2 

tions  for  Ai{k^^,)  and  in  (44)  for  1  corresponding  to  the  dark  side 

of  the  caustic  which  lies  exterior  to  the  caustic  transition  region;  these  asymptotic 
approximations  are: 


Ai{k^,) 


>>  1 


Ai'{klC)^-^{klc)\e- 

2  y'tt 


» 1 


It  is  important  to  observe  that  even  though  the  values  associated  with  both  complex 

stationary  points  Q'  and  Ql,  respectively,  are  present  in  the  expressions  for 

in  (45),  the  use  of  (47)  and  (48)  in  (44)  yields  an  asymptotic  limit  containing  only 
the  field  reflected  from  the  complex  stationary  point  Q\  (and  not 
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Ill  An  Asymptotic  Solution  for  the  Field  Associ¬ 
ated  with  a  Three-Dimensional  Smooth  Caus¬ 


tic: 

The  development  presented  in  Section  2  for  the  2-D  case  is  repeated  here  for  the 
corresponding  3-D  situation.  The  excitation,  as  in  the  2-D  case,  is  assumed  to  be 
produced  by  a  localized  source  such  that  the  field  incident  on  the  3-D  reflecting 
boundary  can  be  represented  ray  optically.  Again,  both  the  non-uniform  and  the 
uniform  asymptotic  approximations  of  the  PO  integral  for  the  scattered  field  in 
.3-D  are  presented  separately  in  Parts  A  and  B  of  this  section,  respectively.  As 
before,  the  non-uniform  asymptotics  yield  the  pair  of  GO  reflected  ray  fields  on 
the  lit  side  of  the  smooth  caustics  of  reflected  rays  that  are  associated  with  the 
two  real  stationary  points  correspond  to  points  of  reflection  at  Qa  and  Qi,  on  the 
3-D  reflecting  boundary.  These  two  reflected  GO  rays  coalesce  at  the  caustic  where 
their  non-uniform  GO  field  becomes  unbounded.  The  non-uniform  field  on  the  dark 
side  of  the  caustic  is  associated  with  only  one  of  the  two  complex  stationary  points 
which  yields  an  exponentially  decaying  field  behavior  (the  other  does  not  satisfy 
the  radiation  condition);  this  complex  stationary  point  is  denoted  by  Ql  on  the 
complex  extension  of  the  reflecting  boundary.  Once  again,  it  is  noted  that  the  non- 
uniform  analysis  is  not  valid  at  and  near  the  smooth  caustics  of  reflected  rays;  one 
must  instead  employ  the  results  based  on  the  uniform  analysis  within  these  caustic 
regions.  The  non-uniform  result  for  the  dark  side  of  the  caustic  is  also  unbounded  at 
the  caustic.  Of  course,  the  uniform  results  reduce  to  the  non-uniform  results  outside 
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the  caustic  transition  regions.  As  done  for  the  2-D  case,  the  non-uniform  analysis 
serves  to  introduce  the  notation  for  the  3-D  case  which  will  also  be  employed  in  the 
uniform  asymptotic  analysis. 

In  the  3-D  case,  there  can  in  general  be  two  smooth  caustics  of  reflected  rays; 
the  present  analysis  is  valid  across  each  caustic  provided  that  these  two  caustic 
surfaces  do  not  intersect.  For  the  sake  of  convenience  of  analysis,  it  is  assumed 
that  the  3-D  reflecting  boundary  is  rotationally  symmetric  (about  the  z-axis  (say)), 
and  that  the  observation  point  lies  in  the  plane  defined  by  the  point  source  (which 
illuminates  the  boundary)  and  the  axis  of  revolution  (or  z-axis)  of  the  rotationally 
symmetric  boundary.  Under  this  assumption,  the  plane  of  incidence  (or  reflection) 
defined  by  the  incident  ray  direction  and  the  normal  to  the  surface  both  of  which 
are  evaluated  at  the  real  point  of  reflection  on  the  boundary  coincides  with  one  of 
the  principal  planes  of  that  boundary.  While  this  situation  arises  only  in  special 
cases,  it  is  chosen  because  it  simplifies  the  analysis.  Furthermore  the  asymptotic 
results  obtained  in  this  special  case  can  be  employed  directly  to  deal  also  with  the 
general  case  in  which  the  boundaries  are  not  necessarily  rotationally  symmetric, 
and  in  which  the  planes  of  incidence  do  not  necessarily  coincide  with  any  of  the 
principal  planes  of  the  surface  at  real  points  of  specular  reflection  on  the  boundary. 

III-A  Non-Uniform  3-D  Analysis: 

The  electric  field  E*(P)  which  is  scattered  by  a  3-D  boundary  illuminated  by  a 
ray  optical  field  (chosen  here  to  be  that  due  to  a  point  source  at  P'),  and  which  is 
observed  at  a  point  P  exterior  to  the  boundary  can  be  expressed  as  usual  in  terms 
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of  the  induced  current  J  on  the  boundary  S  by: 


jj\kx  k< 


in  which  Q'  is  any  point  on  the  boundary  S  and  R  is  the  vector  from  Q'  to  P  (or 
R  =  Q'P  )  as  before.  Also,  P  cannot  be  too  close  to  5  for  (49)  to  be  valid.  The 
PO  approximation  to  (49)  is  employed  via  the  usual  GO  approximation  J{Q')  — 
2h  X  H'{Q'),  (on  the  lit  portion  of  S  and  zero  otherwise)  given  in  (2),  where  h  is 
the  unit  outward  normal  vector  to  5  at  Q'  and  the  incident  ray  optical  magnetic 
field  H'{Q')  due  to  a  point  source  at  P'  is  given  by 

H'{Q')  -  X  A  .  ;  s’  .  A  =  0  (50) 

s* 

As  usual,  H*{Q')  w  Z~^3'  x  E'{Q'),  in  which  the  incident  electric  field  E'{Q')  = 
A  *  .  It  is  noted  that  s’  =  P'Q'  as  in  the  2-D  case,  and  s*  =  . 

Since  one  is  primarily  interested  in  the  fields  specularly  reflected  from  the  surface, 
it  will  be  assumed  that  the  boundaries  of  the  PO  integration  (resulting  from  the 
shadow  boundary  on  the  surface)  are  far  from  the  stationary  points  in  the  expo¬ 
nential  phase  term  of  the  integrand;  furthermore,  these  spurious  contributions  to 
the  PO  integral  arising  from  the  shadow  boundary  will  be  ignored.  For  the  sake  of 
convenience,  the  region  of  integration  will  still  be  denoted  by  S.  Thus,  under  the 
PO  approximation,  (49)  becomes: 


E’{P) 


k  X  R  X  {h  X  i'  X  A)- 


-jk{$'+R) 


Let  ds'  =  di'  •  dTj,  where  is  the  coordinate  curve  corresponding  to  the  gener¬ 
ator  of  the  3-D  surface  of  revolution  S  and  Tj  is  the  azimuthal  coordinate  which  is 
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orthogonal  to  As  in  the  2-D  case,  the  phase  term  exhibits  real  station¬ 

ary  phase  points  at  Qa  and  on  S  '\{  P  —  Pi  on.  the  lit  side  of  the  smooth  caustic; 
likewise  it  exhibits  complex  stationary  points  at  Q'  and  Q%  (which  are  complex 
conjugates  of  each  other)  if  P  =  Ps  oii  the  dark  side  of  the  caustic.  The  stationary 
phase  condition  is  given  by: 

t\  ■  V(s*  +  P)  =  0  ;  Pj  •  V(3’  +  P)  =  0  .  (52) 

Both  the  above  conditions  (namely  (52a)  and  (52b))  must  be  satisfied  simultane¬ 
ously  at  the  stationary  points.  If  the  source  point  and  the  axis  of  revolution  of  S 
define  the  plane  and  if  the  observation  point  lies  in  this  plane  then  (52a; 

52b)  can  be  shown  to  yield  stationary  phase  points  for  which 

^2  ~  ^2«  stationary  phase  points)  (53) 

along  with  i\  =  for  Qa  and  i\  =  at  Qj,,  if  P  =  P^.  Thus,  the  coordinates  of 
Qa  are  (^lal^i*)  of  respectively.  Likewise,  when  P  =  Ps,  the 

coordinates  of  Q^  are  of  Qb  quantities  PiaJ^isj^ia 

and  have  the  same  meaning  as  in  (6a;  6b)  and  (7)  for  the  2-D  case. 

At  the  real  stationary  points  Qa  and  Qj,  on  S  for  P  =  Pi  the  realtionships 
given  in  (8)-(ll)  for  the  2-D  case  also  hold  for  the  present  3-D  situation  because  the 
condition  "=  ^2»  stationary  phase  points  describes  a  plane  which  reduces 

the  3-D  problem  effectively  into  a  2-D  problem.  The  reflected  rays  in  Figures  1(a) 
and  1(b)  for  the  2-D  case  also  exist  in  the  plane  Pj  =  ^2,  for  the  3-D  case  except  that 
the  boundary  in  those  figures  must  be  modified  to  become  rotationally  symmetric 
for  the  present  3-D  analysis. 
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In  the  neighborhood  of  the  stationary  points  and  Q5  when  P  =  Pl,  one  may 
employ  the  following  local  approximation  for  the  phase  term  /c(s*  -|-  R): 

,<  +  R^  +<,i+  -  (Im)’ 

for  Q'  near  Qa,b  ■  (54) 

At  high  frequencies  (or  large  k  ),  the  dominant  contribution  to  (51)  comes  from 
the  stationary  phase  points  and  Qb  on  5  when  P  =  Pt  \  therefore,  (51)  is 
approximated  via  (54)  to  obtain: 

^  ^  {[«  X  X  (n  X  i*  X  A{Q'))]  .  ^ 

y  |{^  X  ^  X  (n  X  i*  X  A(Q'))]~ 


}]  (55) 


in  which  Gi{Qa,b)  is  given  as  before  in  (14)  together  with  (15).  Likewise,  G2{Qa,b) 


is  defined  as: 


G2iQa,b)  —  ~  t  r  ’  Qo>> 
^a,b  P2a,2b 


It  can  be  shown  that 


1  2  cos  0' 


p2a,2b  ^a,b  PtiQa.b) 

The  domains  of  integration  si  and  S2  in  (-55)  include  the  disjoint  neighborhoods  of 
the  points  and  Qf,  on  S.  In  (57),  Pi(Qa.b)  denotes  the  radius  of  curvature  of  S 
along  at  Qa,b  ;  whereas  Pg{Qa,b)  which  occurs  in  (14)  via  (15)  is  the  radius  of 
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curvature  of  S  along  at  Qa,b  as  indicated  previously.  Additional  approximations 
valid  near  Qa  and  Qi,  in  the  integrals  over  Si  and  32  of  (55)  yield 


-  1  h  i  ~ 

E‘{Pl)  ^  <  [fia  X  -Ra  X  (h„  X  j;  X  A(Q„))] - — — 

Ztt  [  a; 

r  df,  •  r  dC  e-J^* <‘5- 

7-00  ^  7-00 

+  [kb  X  RbX  (hfc  X  i;  X  i(Qfc))]^ — 

.  r  di\  •  r  de[  .(58) 

J-oo  J —00  ) 

Following  essentially  the  same  steps  which  lead  one  from  (15)  to  (16)  for  the  2-D 
case,  it  can  be  easily  seen  that  (56)  reduces  in  the  non-uniform  stationary  phase 
evaluation  to: 

E-(Pl)  «  ^  hfi,  X  fl,  X  (ft,  X  i;  X  4(Q„))] - p;— 


kG  liQ  a)  cos^Oi 


,-jfSffn(G,  ((?.))  ,  /  g->?50n[G,(Q.)] 

\lkG2{Qa) 


+  [Ri,  X  «{,  X  (ttj,  X  i;  X  i((?{,))]^ — 


/  j  ’Sflnl<3l(C?»)l 

1  ^TT  -j;S|,n[G2«?t)ll 

^  kG,{Qb}cos^  01 

V  kG2{Qb)  J 

Making  use  of  the  previously  mentioned  relationships  in  (14),  (17),  (21),  (24)  and 
(56)  simplifies  (59)  to  the  more  familiar  GO  reflected  field  from  as: 


E'(n)  ~  ■ ««?,)./  7^ 

V  Pla  +  -^o  V  P2a  + 

+E'(Qb)  •  k{Qb)J 

V  P\b  +  V  P2b  +  ^b 
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It  is  noted  that  E‘{Pi)  in  (60)  is  the  sum  of  the  GO  fields  refiected  from  the  points 
Qa  and  Qt,  on  5  to  Pi  on  the  lit  side  of  the  caustic  defined  by  where 

.  Also,  the  relationships  in  (22)  and  (23)  are  true  in  the  result  of  (60); 

furthermore, 

if  S  is  a  surface  of  revolution  because  2(>  >  G2(Qa,b)  >  0  in  that  case  since 

the  axis  of  revolution  is  now  the  other  caustic  of  the  reflected  rays  (i.e.,  the  second 
caustic  surface  degenerates  to  a  line). 

The  result  in  (60)  fails  at  and  near  the  caustics  of  reflected  rays.  Although 
(60)  has  been  developed  for  the  special  case  of  a  surface  of  revolution  with  the 
observation  point  Pi  being  located  in  the  plane  defined  by  the  source  point  and  the 
axis  of  revolution,  it  is  also  valid  for  the  GO  reflection  of  an  arbitrary  ray  optica] 
field  from  an  arbitrary  surface,  except  that  and  longer 

be  obtained  from  the  simple  expressions  in  (15)  and  (57),  respectively,  but  instead 
they  must  now  be  obtained  from  the  more  complicated  expressions  in  [16,17,18]  for 


the  latter  more  general  situation. 

A  similar  stationary  phase  analysis  at  a  point  P  =  Ps  on  the  shadow  side  of  the 
caustic  defined  by  sj  =  with  p\^,  =  -  yields  the  following  expression  for 
the  scattered  field  E*{Ps)' 


E‘(Ps)-E'{Ql) 


(62a) 


where, 

E'iQl)  -  A((?^)"  ,  ‘ 
•’I 


(62b) 
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The  expressions  in  (62a;  62b)  are  analogous  to  those  in  (26a;  26b)  for  the  2-D  case. 
As  mentioned  earlier,  the  field  in  (62a)  is  specularly  reflected  from  the  point 
which  lies  on  the  complex  extension  of  the  original  boundary;  the  contribution  from 
the  other  (complex  conjugate)  stationary  phase  point  at  Q®  must  again  be  discarded 
on  physical  grounds  as  it  violates  the  radiation  condition.  The  contribution  in  (62a) 
exhibits  an  exponentially  decaying  field  at  F5.  As  in  the  2-D  case,  all  the  quantities 
associated  with  the  superscript  c  must  be  evaluated  via  an  analytic  continuation 
procedure.  The  non-uniform  stationary  phase  approximation  in  (62a)  is  not  valid 
at  the  caustic  associated  with  pjj,  ;  it  is  unbounded  there.  A  uniform  asymptotic 
approximation  which  remains  valid  at  and  near  this  caustic  is  presented  below. 

III-B  Uniform  3-D  Analysis 

One  begins  the  uniform  asymptotic  analysis  by  noting  that  the  non-uniform  results 
in  (60)  and  (62),  which  become  unbounded  at  sj  =  with  p\^  =  —  |pi^|,  are  also 
valid  for  the  reflection  of  an  arbitrary  ray  optical  incident  field  from  an  arbitrary 
.3-D  surface.  Furthermore,  the  only  difference  between  the  non-uniform  results  in 
(60)  and  (62),  for  the  3-D  case  and  the  corresponding  results  in  (25)  and  (26)  for  the 
2-D  case  are  the  additional  reflected  ray  divergence  factors  ./ and 
which  are  present  in  (60)  and  (62)  respectively  for  the  3-D  case.  The  3-D  case 
reduces  to  the  2-D  case  when  26  \Pv>\  — >  00  so  that  the  spreading  of 

the  reflected  rays  is  now  restricted  only  to  one  plane  for  the  2-D  case,  and  the  other 
spreading  (or  ray  divergence)  factors  alluded  to  above  (containing  P2a.2i>  PVb) 
become  unity  for  the  2-D  case.  Consequently,  the  uniform  GO  result  for  the  3-D 
case  which  reduces  to  (60)  and  (62)  on  the  lit  and  dark  sides,  respectively,  outside 
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the  caustic  transition  regions  can  be  obtained  directly  from  (33)  and  (44)  together 
with  (37)  and  (45),  respectively,  by  conjecturing  that  for  the  lit  side,  E‘{Pl)  is  given 

by 

E‘{Pl)  ~  2n  Qt  k^"^  Ai'i-kUt)]  (63) 


with  6(  and  Q  (where  >  0  )  being  the  same  as  in  (28)  and  (29).  Also,  Pi  and 
Q(  are  given  by  a  simple  modification  of  (37)  to  include  the  additional  incident  and 
reflected  ray  spatial  spreading  (or  divergence)  factors;  thus. 


f  A  1 

U-J 

2  Vtt 

i(C?J  •«((?„)  / 


P2a 


.s!. 


P'la  +  \  P2a  + 


A{Qk)-R{Q6)  I  P\, 


P2h 


+  4  ^P2b  +  4 


(64) 


Likewise,  E‘(Ps)  for  the  shadow  side  is  given  by 


E’iPs)  -  2irj  e-^  ''"•(P.fc-i  At(A:^C.)  +  k^i  Ai'i-k^C,)] 


1  -  i 


(65) 


with  S,  and  (where  >  0  )  being  the  same  as  in  (41)  and  (42).  In  addition,  the 
P,  and  Q,  are  given  by 


P. 

Q. 


MQDRiQl) 


PVa 


PZ 


PVa  +  V  PVa  + 


A{Ql).R{Ql)  I  p\l 


PVb 


PVb  +  V  PVb  +  ^6 


(66) 


Thus,  the  Pi  s  and  Qis  above  in  (64)  and  (66)  are  obtained  by  simply  including 

the  transverse  ray  divergence  factors  ./  and  —k—  ./  into  the 

V’L*  V  V*«%  V 

2-D  results  given  in  (37)  and  (45);  here  is  the  additional  transverse  spreading 
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or  divergence  factor  for  the  incident  field  associated  with  a  point  source,  whereas 
is  the  additional  spreading  factor  for  the  reflected  ray  in  the  3-D  case. 
The  results  in  (63)  and  (65)  are  also  valid  for  the  general  3-D  case  as  are  the  non- 
uniform  results  in  (60)  and  (62),  with  the  understanding  that  P2a,2b  ) 

i P'ia, ib'i  P 2a, 2b)  must  be  found  from  the  more  general  expressions  in  [16,17,18]  that  are 
valid  even  when  the  plane  of  incidence  does  not  coincide  with  any  of  the  principal 
planes  of  the  surface  at  the  point  of  reflection.  Furthermore,  the  illumination  (or 
the  incident  field)  may  be  an  arbitrary  ray  optical  field  (rather  than  a  special  point 
source  type  illumination).  For  the  point  source  or  spherical  wave  illumination, 

arbitrary  ray  optical  incident  field,  the  terms 
in  (64)  and  (66)  must  be  replaced  by  E\fiQ{Q*)  where  E\ng(Q')  is 

the  incident  ray  optical  field  whose  phase  at  Q'  is  given  by  exp(— jfcs^;js(Q')).  It  is 
easily  verified  that  the  results  in  (63)  and  (65)  which  are  valid  at  the  caustic  defined 
with  pjj,  =  —  Ipjjl  reduce  uniformly  to  the  results  in  (60)  and  (62), 
respectively,  outside  the  caustic  transition  regions  where  0  >>  1  and  >>  1;  these 
limits  leading  to  (60)  and  (62)  can  be  obtained  via  the  asymptotic  approximations 
for  the  Airy  functions  described  previously  in  (38;39)  and  (47;48). 

The  uniform  results  conjectured  in  (63;64)  for  P  =  Pi  could  be  justified  by 
starting  with  the  PO  integral  for  E*{P)  in  (51)  by  introducing  a  local  approximation 
analogous  to  (27)  for  the  phase  term  in  the  integrand  of  (51)  which  unlike 

(54)  is  valid  even  when  and  Qi,  are  close  together;  namely 

s'  +  R^  ^  P(([)  0  HO  4  0  +  \  G2{nO)i^2  -  Of  (6T) 

with  and  Q  as  in  (28)  and  (29)  for  the  special  case  when  the  observation  point  at 
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P  =  Pi  lies  in  the  plane  defined  by  the  source  point  and  the  axis  of  revolution  as 
was  assumed  in  the  previous  non-uniform  stationary  phase  analysis  for  (51)  which 
led  to  (60).  The  G2{t{([))  is  evaluated  at  =  ^2,  and  has  the  general  form  as  in 
(56);  namely 


1  1 


(68) 


with  a’'(f'i„ifc)  =  and  =  P2a.2b  >  which  and  i\b,^2,  are  the 

coordinates  of  the  stationary  phase  points  Qa  and  Qi,  on  the  line  as  indicated 

below  (53).  Incorporating  (67)  into  (51)  yields: 


E‘(P)^  e--'**'  r  dt 

J  —  OO 

r  (69) 

J -<5o  V  27r 

with 

(?((,<',)  s  X  R  X  (n  X  .i'  X  ^  ^  (70) 

where  f  =  t{i\)  as  noted  earlier.  Also,  the  region  of  integration  has  been  extended 
to  infinity  via  the  usual  stationary  phase  arguments.  Evaluating  the  inner  integral 
over  first  using  the  stationary  phase  approximation  yields 

r  it  I  (71) 

■'-»  I  lAW)  I 

In  obtaining  (71),  the  quantity  )  in  (69)  is  approximated  by  its  value  at  the 

stationary  point  when  evaluating  the  inner  integral  on  (.'2  in  (69).  The 

remaining  integral  in  (71)  on  /(f'j)  can  be  evaluated  exactly  as  in  (30)-(33)  for  the 
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2-D  case  to  obtain  the  form  of  the  solution  which  is  the  same  as  in  (63)  with  Pt  and 
Q(  given  by 


-1* 


47  k  4(Q,).fl(QJcos^; 

2  V  TT  4^;  sjGAQa)  cos;  ^G2{Qa) 


A(Q,)  ■  R{Q,)  cos  91 


(72) 


\/Gi{Qb)  cos;  ^G2{Qb) 

T 

in  which  G'j''(Qa,i))  =  ,  and  ^'^'(Qa.fc)  is  the  same  as  G7^(Qa,b)  with  the 

subscript  1  replaced  by  2.  It  is  then  easily  verified  in  the  special  case  of  a  surface 
revolution,  with  Pi  being  located  in  the  plane  containing  the  axis  of  revolution  of 
the  surface  and  the  source  point,  that  (72)  is  indeed  identical  to  the  one  conjectured 
in  (64)  using  physical  arguments. 

The  result  conjectured  in  (65;66)  for  the  dark  side  of  a  caustic  can  also  be 
justified  similarly  for  the  special  case  of  the  surface  of  revolution  when  Ps  is  located 
in  the  plane  containing  the  axis  of  revolution  and  the  source  point.  For  the  latter 
case,  the  6(  and  —4/  in  (67)  must  be  replaced  by  6,  and  +4,  w'here  4s  >  0,  with  6, 
and  4s  is  as  defined  in  (41)  and  (42),  respectively. 


IV  Numerical  Results 

The  2-D  uniform  asymptotic  solution  given  in  (33)  and  (37)  has  been  applied  here 
to  study  the  backscattering  from  boundaries  which  contain  points  of  inflection. 
One  such  application  is  shown  in  Figure  2  where  a  closed  concave-convex  shaped 
geonretry  with  an  edge  is  illuminated  by  a  plane  wave,  and  it  is  of  interest  to  find 
the  far  zone  backscattered  field.  In  addition  to  the  field  of  the  rays  specularly 
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reflected  from  this  boundary,  the  edge  diffracted  ray  fleld  is  also  included  in  this 
backsca*ter  calculation;  however,  the  effect  of  the  rays  that  creep  around  the  convex 
portion  of  the  surface  after  being  excited  directly  at  grazing  incidence  or  by  edge 
diffraction  is  ignored  in  this  work.  Also  neglected  are  the  effects  connected  with 
the  launching  of  whispering  gallery  type  modes  on  the  concave  boundary  by  the 
mechanism  of  edge  diffraction,  and  their  subsequent  transformation  into  creeping 
wave  or  surface  ray  modes  on  the  convex  position  of  the  boundary  past  the  inflection 
point.  The  latter  effects  which  have  been  ignored  may  not  be  negligible  far  from 
the  reflected  ray  caustic  transition  region;  hence,  the  backscattered  field  has  been 
calculated  using  the  uniform  GO  results  of  (33)  and  (37)  together  with  the  singly 
edge  diffracted  contribution  only  for  20°  <  0  <  80°  in  Figures  3  and  4  for  the  case 
when  the  incident  electric  field  is  perpendicular  (i.e.,  soft  case)  and  parallel  (i.e., 
hard  case)  to  the  plane  of  incidence,  respectively.  Also  shown  in  Figures  3  and  4 
are  comparisons  with  an  independent  moment  method  (MM)  based  calculation;  the 
good  agreement  between  the  uniform  solution  and  the  corresponding  MM  solution 
serves  to  confirm  the  accuracy  of  the  uniform  solution  on  both  the  lit  and  the  dark 
side  of  the  caustic  of  reflected  rays.  It  is  important  to  note  that  the  maximum  in 
the  backscattered  pattern  occurs  near  B  —  .50°  which  is  just  on  the  lit  side  of  the 
caustic  of  reflected  rays.  Thus,  the  effect  of  the  ray  caustic  is  quite  significant  in 
producing  a  strong  backscatter  around  this  aspect  angle  of  ^  =  50°.  It  is  further 
noted  that  the  region  corresponding  to  B  <  .50°  is  on  the  dark  side  of  the  caustic. 

Additional  calculations  via  the  present  uniform  GO  analysis  have  been  per¬ 
formed  for  the  far  zone  backscattering  by  a  2-D  smoothly  indented  cavity  of  Figure 
5,  which  is  illuminated  by  an  electromagnetic  plane  wave.  The  results  of  these  cal- 
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(  Illations  are  presented  in  Figures  7-11.  Actually  there  are  three  parts  to  Figures 
7-9;  in  part  (a)  of  these  figures  the  result  for  the  backscattering  from  the  geometry 
of  Figure  5  are  presented,  whereas,  in  parts  (b)  and  (c)  of  each  of  these  figures 
the  backscattering  from  a  2-D  rectangular  cavity  of  the  same  length  and  depth 
are  presented  for  the  sake  of  comparison  [19].  The  2-D  recta  igidar  cavity  geome¬ 
try  is  shown  in  Figure  6.  The  results  in  parts  (b)  and  (c)  of  Figures  7-9  for  the 
backscattering  by  a  rectangular  cavity  pertain  to  the  cases  when  the  incident  elec¬ 
tric  field  is  parallel  and  perpendicular  to  the  plane  of  incidence,  respectively.  The 
results  in  part  (a)  of  Figures  7-9  are  the  same  for  both  of  these  types  of  incident 
field  polarization  over  the  range  of  aspects  |^|  <  60°  which  have  been  considered. 
The  latter  is  true  because  the  polarization  sensitive  creeping  or  surface  ray  modes 
launched  on  the  convex  portion  of  the  reflecting  boundary  outside  the  smoothly 
indented  cavity  portion  do  not  contribute  in  the  range  j^|  <  60°,  and  because  the 
uniform  GO  reflection  contribution  from  the  smooth  indentation  gives  the  same  dB 
magnitude  levels  for  the  two  different  incident  polarizations  (only  the  phase  of  the 
total  uniform  GO  field  differs  by  a  sign  for  the  two  types  of  incident  polarization). 
Figures  7-9  pertain  to  cavities  of  the  same  length;  only  the  cavity  depth  changes 
in  each  of  these  figures.  Figures  10  and  11  show  the  change  in  the  backscattering 
from  the  smoothly  indented  cavity  of  Figure  5  due  to  changing  the  cavity  length 
while  keeping  the  cavity  depth  fixed.  It  is  seen  from  Figures  7-11  that  the  pertur¬ 
bation  in  the  backscattering  due  to  the  presence  of  a  cavity  is  limited  to  a  smaller 
range  of  aspects  for  the  smoothly  indented  cavity  as  compared  to  that  for  the  sharj) 
edged  rectangular  cavity.  This  is  to  be  expected  because  the  effect  of  the  smoothly 
indented  cavity  is  significant  only  within  the  two  reflected  ray  caustic  transition 
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regions.  Nevertheless,  it  is  also  observed  that  the  backscattering  by  the  smoothly 
indented  cavities  can  be  approximately  of  the  same  order  as  the  rectangular  cavity 
for  aspects  that  lie  between  the  reflected  ray  caustics.  It  is  noted  that  there  are 
two  reflected  ray  caustics  in  the  far  zone  for  the  problem  in  Figure  5  (there  is  one 
caustic  associated  with  each  of  the  two  points  of  inflection  on  the  smooth  cavity). 
These  caustics  are  non-intersecting  in  the  far  zone;  however,  they  can  intersect  to 
form  higher  order  caustics  in  the  near  zone.  The  present  uniform  GO  analysis  can 
be  employed  to  And  the  scattering  from  the  smoothly  indented  cavity  of  Figure  5 
only  for  cavity  lengths  which  are  large  enough  so  that  the  two  points  of  inflection 
and  their  associated  reflected  ray  caustic  transition  regions  do  not  come  close  to¬ 
gether.  Finally,  the  present  PO  based  uniform  analysis  also  cannot  be  used  near 
grazing  angles  of  incidence  on  the  reflecting  boundaries  with  points  of  inflection.  It 
is  proposed  as  a  part  of  future  effort  to  remove  the  latter  restriction.  In  conclusion, 
the  present  uniform  GO  solution  is  shown  to  be  accurate  and  relatively  simple  to 
use  within  its  domain  of  validity;  the  conditions  under  which  the  uniform  GO  re¬ 
sults  remain  valid  are  essentially  no  different  than  the  conditions  required  for  the 
conventional  GO  analysis  to  be  valid. 
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Appendix 


The  Parameters  Used  in  the  Uniform  GO  Calcu¬ 
lation 

In  this  appendix,  one  defines  the  parameters  which  have  been  used  in  the  uniform 
GO  solution  for  the  lit  side  of  the  caustic.  These  parameters  are  then  analytically 
continued  to  complex  coordinate  space  in  order  to  calculate  the  field  on  the  shadow 
side  of  the  caustic.  These  parameters  are  presented  for  the  general  3-D  case;  however 
one  can  easily  obtain  the  2-D  results  by  requiring  2  =  0  in  the  3-D  results. 


Parameters  for  the  Lit  Side  of  the  Caustic 

Let  the  coordinates  of  the  source  point  P'  and  the  observation  point  Pi  be  denoted 
by  {x',  y',z')  and  (xi,yi,  zi),  respectively.  For  the  sake  of  being  specific,  one  may 
choose  the  real  point  of  reflection  Qb  on  the  actual  reflecting  boundary;  let  the 
coordinates  of  this  point  Qb  on  the  surface  be  (xb,yb,  Zb)-  Then  one  can  define  the 
ray  path  lengths  $1,  sj,  the  direction  of  the  incident  ray  5^,  the  direction  of  the 
reflected  ray  ij,  the  surface  normal  rib  at  Qb,  the  scalar  product  between  the  surface 
normal  and  the  reflected  ray  direction  cos  01,  the  radius  of  curvature  of  the  surface 
Pg  at  Qb  and  the  vectors  e^,  c[|,  and  ejj  associated  with  the  incident  and  reflected 
ray  for  the  point  of  Qb  as  the  following: 


4  =  \J{^b  -  x')^  +  {yb  -  y'Y  +  (2b  -  z'Y  (73) 
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cos(2^b)  = 
cos  = 

hb  = 


yj{xi  -  Xbf  +  (Vl  -  Vbf  +  {ZL  ~  Zbf 
i  [(i6  -  x')x  4-  ivb  -  y')y  +  [zb  -  z')z] 

4 

—  [{xL  -  Xb)x  +  (y^  -  yb)y  +  {zl  -  Zb)z] 

4 


-■Si 


lcos{20l)  +  1 
2  . 


H  - 


2  cos 


"^  =  in4i 

ej|  =  ex  X  i|, 
e[|  =  ex  X 

Pg(Qb)  =  Pg(Xb,yb,Zb) 


(74) 

(75) 

(76) 

(77) 

(78) 

(79) 


(80) 


Parameters  for  the  Shadow  Side  of  the  Caustic 


When  the  observation  point  is  on  the  shadow  side  of  the  caustic  the  point  of  reflec¬ 
tion  is  no  longer  real  (as  is  Qb)i  is  now  replaced  by  a  complex  point  of  reflection 
Qb  as  discussed  in  Sections  II  and  III  of  the  report.  The  corresponding  parameters 
for  the  shadow  side  are  denoted  by  the  same  symbols  as  in  (74)-(81)  for  the  lit 
side  except  that  they  now  contain  a  superscript  c  to  denote  the  fact  that  they  are 
associated  with  the  complex  point  Qj.  Let  the  coordinates  of  the  observation  point 
Ps  on  the  dark  side  be  (x,,y,,2,).  Let  the  coordinate  of  the  complex  reflection 
point  Qb  on  the  complex  (or  analytically  continued  part)  of  the  extended  surface 
be  (Xb,yb,  zp.  One  can  then  redefine  the  above  parameters  for  the  lit  side  so  that 
they  will  be  valid  for  the  shadow  side  by  simply  replacing  the  coordinates  of  Qb  by 
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those  of  Q\.  Thus,  one  arrives  at  the  following  results: 


cos(2^n  = 

cos  = 

hi  = 


+(y^^_y/)2+(^c  _  ^,)2 

\J{XL  -  xl)^  +  (t/i  -  yt)^  +  {zl  -  zt)^ 

[(®6  -  a;')*  +  ivl  -  y')y  +  {4  -  ^')4 

-4  “  4)^  +  (y^  -  yt)y  +  -  zi)z] 

^6 


(81) 

(82) 

(83) 

(84) 

(85) 

(86) 

(87) 

(88) 
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A  CONCAVE  BOUNDARY 


(b)  REFLECTION  FROM  A 
CONCAVE  -  CONVEX 
BOUNDARY 

1:  Scattering  by  smooth  boundaries  which  produce  smooth  caustics  of 
rays. 


re- 


0 


X 


The  boundary  is  defined  by  the  equation 

y=  i>iyf(cos(^+l) 

with  A  =  0.3  m  being  the  maximum  length  in  the  x  direction 
and  the  frequency  at  which  the  calculations  of  backscatter  are  made  is  3  GHz. 

Figure  2:  Scattering  by  a  boundary  with  a  point  of  inflection. 
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Figure  3:  Backscattering  from  the  geometry  in  Figure  2  when  the  electric  field  is 
polarized  perpendicular  to  the  plane  of  incidence. 
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Figure  4:  Backscattering  from  the  geometry  in  Figui 
polarized  parallel  to  the  plane  of  incidence. 
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Figure  5:  Smoothly  indented  cavity  configuration. 
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ure  6:  Rectangular  cavity  configuration. 


cavity  (Figure  5)  of  L  =  10.2A  and  D  =  l.A. 

Figure  7:  Backscattering  by  a  smoothly  indented  cavity  (of  Figure  5)  as  well  as  by 
a  rectangular  cavity  (of  Figure  6). 
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ANGLE  IN  DEGREES 

(b)  Backscattering  by  a  rectangular  cavity  (Figure  6) 
of  £  =  10.2A  and  D  =  l.X  when  the  electric  field  is 
parallel  to  the  plane  of  incidence. 

Figure  7:  (continued) 
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ANGLE  IN  DEGREES 
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(c)  Backscattering  by  a  rectangular  cavity  (Figure  6) 
oi  L  =  10. 2A  and  £1  =  l.A  when  the  electric  field  is 
perpendicular  to  the  plane  of  incidence. 

Figure  7:  (continued) 
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(a)  Backseat tering  by  a  smoothly  indented 


cavity  (Figure  5)  o{  L  —  10.2A  and  D  =  I.IA. 


Figure  8:  Backscattering  by  a  smoothly  indented  cavity  (of  Figure  5)  as  well  as  by 
a  rectangular  cavity  (of  Figure  6). 
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ANGLE  IN  DEGREES 

(b)  Backscattering  by  a  rectangular  cavity  (Figure  6) 
of  £  =  10. 2A  and  D  =  1.1  A  when  the  electric  field  is 
parallel  to  the  plane  of  incidence. 

Figure  8:  (continued) 
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ANGLE  IN  DEGREES 

(c)  Backscattering  by  a  rectangular  cavity  (Figure  6) 
of  £  =  10. 2A  and  D  =  I.IA  when  the  electric  field  is 
perpendicular  to  the  plane  of  incidence. 

Figure  8:  (continued) 


•57 


MAGNITUDE  IN  DB  PHASE  IN  DEGREES 
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(a)  Backscattering  by  a  smoothly  indented 


cavity  (Figure  5)  of  L  =  10.2A  and  D  —  1.25A. 


Figure  9:  Backscattering  by  a  smoothly  indented  cavity  (of  Figure  5)  as  well  as  by 
a  rectangular  cavity  (of  Figure  6). 
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ANGLE  IN  DEGREES 

(b)  Backscattering  by  a  rectangular  cavity  (Figure  6) 
of  X/  =  !0.2A  and  D  =  1.25 A  when  the  electric  field  is 
parallel  to  the  plane  of  incidence. 

Figure  9:  (continued) 
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(c)  Backscattering  by  a  rectangular  cavity  (Figure  6) 
oi  L  =  10. 2A  and  D  =  1.25A  when  the  electric  field  is 
perpendicular  to  the  plane  of  incidence. 

Figure  9:  (continued) 
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Figure  10:  Backscattering  by  a  smoothly  indented  cavity  of  Figure  5  with  L  =  10. A 
and  D  =  l.A. 
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Figure  11:  Backscattering  by  a  smoothly  indented  cavity  of  Figure  5  with  L  =  20. A 


and  D  =  l.A. 
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